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Abstract 

Kinetic energies of light fragments {A < 10) from the decay of target specta- 
tors in Au + 19"^ Au collisions at 1000 MeV per nucleon have been measured 
with high-resolution telescopes at backward angles. Except for protons and 
apart from the observed evaporation components, the kinetic-energy spec- 
tra exhibit slope temperatures of about 17 MeV, independent of the particle 
species, but not corresponding to the thermal or chemical degrees of freedom 
at breakup. It is suggested that these slope temperatures may reflect the 
intrinsic Fermi motion and thus the bulk density of the spectator system at 
the instant of becoming unstable. 



PACS numbers: 25.70.Pq, 21.65. +f, 25.70.Mn 
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Statistical multifragmentation models have been found to be remarkably successful in 
describing the partition space populated in nuclear reactions leading to multi-fragment emis- 
sion This is particularly true for spectator decays following heavy-ion collisions at 
relativistic bombarding energies. Fragment yields and fragment correlations have been re- 
produced with high accuracy [P|-[T3[1, including their dispersions around the mean behavior 
The model temperatures are in good agreement with measured isotope temperatures, 
and the isotopic yields from which these temperatures are derived are well reproduced |TO|,|T3 



The mean kinetic energies of the produced particles and fragments, on the other hand, are 
not universally accounted for if the same set of model parameters is used [^,n)|. Measured 
kinetic energies or, equivalently, the slope temperatures obtained from Maxwell-Boltzmann 
fits to the kinetic energy spectra have been found to considerably exceed the predicted 
values. This holds even though the effects of Coulomb repulsion and sequential decay cause 
the calculated slope temperatures to be higher than the equilibrium temperatures in the 
model description. For the ^^''Au + ^^^Au reaction studied in this paper, the measured 
isotope temperatures and the calculated equilibrium temperatures are 6 to 8 MeV for most 
of the impact parameter range; the calculated slope temperatures for light particles are 
about 10 MeV while the measured slope temperatures exceed 15 MeV ||1CI|| . 

The qualitative difference between the chemical or thermal equilibrium temperatures 
and the kinetic temperatures has been observed for a variety of reactions and numerous 
reasons for it have been presented [Il5|-|l9|]. Preequilibrium or pre-breakup emissions are likely 
explanations for high-energy components in light-particle spectra 0,|l^,^,^ . In some cases, 
the excess kinetic energies have been ascribed to collective flow even though a characteristic 
proportionality to the fragment mass is not evident in the data [|11|,^,^ . Other authors were 
caused by this situation to more generally question the applicability of statistical descriptions 
for spectator fragmentation . An explanation in terms of the Goldhaber model that is 



favored here has been suggested long ago |T5|J26|] but, so far, has not been generally adopted. 

In this Letter, we analyze kinetic-energy spectra for fragment and light-charged-particle 
emission from target spectators following collisions of ^^^Au -|- ^^^Au at 1000 MeV per 
nucleon. We will show that the deduced slope temperatures and the much lower chemical 
breakup temperatures can be consistently understood if the intrinsic constituent motion of 
the decaying fermionic system is taken into account. With this interpretation, the slope 
temperatures reflect the bulk density of the system prior to its disintegration, and thus 
provide information on the mechanism of fragment formation. 

Beams of ^^^Au with incident energy 1000 MeV per nucleon from the heavy-ion syn- 
chrotron SIS were used to bombard ^^^Au targets of 25-mg/cm^ areal thickness. As part 
of a larger experimental setup [|Ty,^,^, three high-resolution telescopes were placed at 
backward angles for detecting the products of the target-spectator decay. Each telescope 
consisted of three Si detectors with thickness 50, 300, and 1000 /im and of a 4-cm long 
CsI(Tl) scintillator with photodiode readout and subtended a solid angle of 7.0 msr. 

The products of the projectile decay were measured with the time-of-fiight wall of the 
ALADIN spectrometer ||^ and Zhound was determined event-by-event. The sorting variable 
Zbound is defined as the sum of the atomic numbers Zi of all projectile fragments with Zi > 2. 
Zbound reflects the variation of the charge of the primary spectator system and is therefore 
correlated with the impact parameter of the reaction. Because of the symmetry of the 
collision system, the mean values of Z^mmd for the target and the projectile spectators within 
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the same event class have been assumed to be identicaL First results of these measurements 
have been presented in Ref. 0], including some of the kinetic-energy data that will be 
discussed in the following. 

Energy spectra of protons, measured at 9iab = 150° and sorted into eight bins of Zhound, 
are shown in Fig. 1 (left panel). They are characterized by a hard component and an 
additional soft component that is most clearly identified at large Zbound- Fits using two 
Maxwellians yield temperature parameters of very different magnitude that both are nearly 
independent of Zbound (Fig. 1, right panels). The slope temperature T of about 5 to 7 
MeV of the soft component is of the same order as the measured breakup temperatures |]TU 



Its intensity, at the same time, is strongly correlated with that of heavier fragmentation 



products ||29| which suggests a, perhaps partial, interpretation as evaporation from highly 
excited residual nuclei. For the ^^^Au + ^^C reaction at the same energy, an equilibrium 
proton component with nearly identical properties was reported by Hauger et al. [p^ . 

Two components have also been identified in the spectra of a particles (Fig. 2) and in 
the neutron data measured with LAND [Q, with low-temperature components typical for 



evaporation. The high-temperature components of protons and neutrons seem to have their 
origin not only in the breakup stage but also in the earlier cascading stages of the collision; 
their slope parameters are rather large and, as established for the neutron case vary 
with the bombarding energy. 

The high-temperature component of a particles decreases from T ^ 11 MeV for small 
Zbound to T ?5i 13 MeV for large Zbound (see below). Its mean value of 15 MeV is close to 
the slope temperatures observed for other species with A > 2 which exhibit Maxwellian 
spectra to a good approximation (Fig. 2). Their slope temperatures were obtained from 
single-component fits with three-parameter functions that included a Coulomb potential Vc- 
The results for the interval 20 < Zbound < 60 are shown in Fig. 3. Apart from the protons 
with T 26 MeV, all spectra exhibit temperatures that are narrowly dispersed around a 
mean value of 17 MeV. 

A mass-invariant temperature, at first sight, seems to provide direct evidence for equi- 
libration of the kinetic degrees of freedom. On the other hand. Coulomb effects should 
contribute in proportion to the fragment charge, recoil effects may be important for heavier 
fragments, and the small but finite motion of the target spectator should introduce a collec- 
tive velocity component. In order to estimate the magnitude of these effects, a model study 
was performed which started from initial configurations generated by randomly placing frag- 
ments and light particles in a spherical volume with random momenta corresponding to a 
given temperature. Experimental charge distributions from Ref. and an average density 
p/Po = 0.3 were used where po is the normal nuclear density. N-body Coulomb trajectory 
calculations were then performed and asymptotic slope temperatures were determined by 
fitting, following the same procedure as with the experimental data. The calculations showed 
that the considered effects are relatively small, perhaps of the order of AT/AA ^ 0.2 MeV, 
and that they cancel each other to a good approximation at the backward angles chosen for 
the measurement. The measured slope parameters oiT 17 MeV (Fig. 3) are thus equal 
to the temperatures to be used in a thermal interpretation of the data. 

The model study also showed that the Coulomb repulsion contributes considerably to 
the final kinetic energies. The resulting mean kinetic energies in the center-of-mass frame of 
the decaying system add up to about 35 to 40 MeV, quite consistent with values obtained 
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in earlier experiments where kinetic energies were derived from the widths of momentum 
distributions of projectile fragments In these experiments, the mean kinetic energies were 
found to be approximately constant over the range of fragment atomic number Z < 20. 

It was suggested many years ago by Goldhaber that the product momenta in fast frag- 
mentation processes may have their origin in the nucleonic Fermi motion within the colliding 



nuclei ||2^. He also pointed out that the resulting behavior is indistinguishable from that 



of a thermalized system with rather high temperature. For the Fermi momentum pp ^ 265 
MeV/c of heavier nuclei the corresponding temperature is T = 15 MeV. This is intriguingly 
close to the measured slope temperatures T = 17 MeV. 

Goldhaber's idea, initially formulated for cold nuclei, has been extended to the case of 



expanded fermionic systems at finite temperature by Bauer |T9[. We have used the numeri 



cal solution reported there by assuming that the temperature of the system is given by the 
measured breakup temperatures TneLi, derived from the yields of '^''^He and ^'^Li isotopes 
ni| , PT|P^ . For the breakup density two alternative values p/po = 1.0 and 0.3 were chosen 



which correspond to two significantly different scenarii for the process of fragment produc- 
tion. The higher value is expected for a fast abrasion stage that will induce fragmentation 
of the residual spectators but not instantly affect the nucleon motion within them. The 
lower value is in the range usually assumed for multi-fragment breakups of expanded sys- 
tems [0,0. Here the development of instabilities, as the system enters the spinodal region, 
will equally cause a rapid fragmentation, so that the nucleonic Fermi motion will contribute 
to the fragment kinetic energies. For a homogeneous system at a lower than normal density 
the Fermi motion is reduced and the effect will be smaller. 

The slope temperatures obtained in this way are represented by the lines shown in Fig. 4. 
The comparison with the data is here restricted to the A < 4 isotopes for which the collected 
statistics are sufficient for studying the Zbound dependence. The recoil factor {As— A) / {Ag — l) 
appearing in the Goldhaber formula [|19| , ^ can be safely ignored as the mass number Ag of 
the spectator system is 100 to 150 on average and still about 50 in the bin of smallest Z^ound 
P^ . Qualitatively, the predicted values are close to those observed. In particular, the rise 
with decreasing Z^ound follows as a consequence of the rising breakup temperatures TneLi- 
Apparently, with assumptions as made in the Goldhaber model, the two different types of 
temperatures are mutually consistent. The role of the slope temperatures, consequently, is 
restricted to describing the fragment motion while TneLi is the more suitable observable for 
representing the temperature of the nuclear environment at the breakup stage. 

Recent calculations with transport models, which incorporate Fermi motion, support 
this interpretation The energies of spectator fragments, as measured in the present 

reaction, are well reproduced, and the coexistence of qualitatively different internal (or 
local) temperatures and fragment slope temperatures has been demonstrated in these QMD 
and BUU studies. On the experimental side, the observed similarity of fragment kinetic 
energy spectra, measured with different projectiles at various incident energies, is expected 
from the Goldhaber model [P, p!5| , ^^ , ^4| -|36[| . The fluctuations corresponding to the higher 
slope temperatures may also explain the observed large widths of the intrinsic kinetic-energy 
distributions in ^^^Au fragmentation at 600 MeV per nucleon that are not easily reproduced 
with statistical multifragmentation models p7| . 

A more quantitative comparison with the data shown in Figs. 3 and 4 will favor the 
prediction for p/pQ = 1.0 corresponding to a fast breakup over that for p/po = 0.3 that is 
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reached after a homogeneous expansion. This would not be unexpected in the present case 
of spectator fragmentation following collisions of heavy ions and might indicate a limited 
resilience of nuclear matter to shape deformations |3^. The cited QMD and BUU calcula- 
tions consistently suggest that the fragments are preformed at an early stage of the collision 
(< 50 fm/c) before the system has expanded to typical breakup densities. 

To draw firm conclusions at this time would certainly be premature. The role of sec- 
ondary decays and possibly other effects would have to be considered, and a moderate col- 
lective velocity component as a result of thermal expansion cannot be completely excluded 
P7| , |5D[] , even though it is not explicitly indicated for the present reaction (Fig. 3 and 0). 
The main purpose of this work is to show that, within the Goldhaber picture of a random 
superposition of the nucleon momenta, the larger slope temperatures are fully consistent 
with a breakup at equilibrium temperatures of 6 to 8 MeV, as they are measured and also 
obtained with statistical multifragmentation models. Because of the sensitivity to the mode 
of fragment formation, via the density at which the Fermi motion has been established, this 
interpretation seems rather attractive and deserves further attention. 

The authors would like to thank J. Aichelin, C. Fuchs, T. Gaitanos and H.H. Wolter for 
fruitful discussions. M.B., J. P., and CS. acknowledge the financial support of the Deutsche 
Forschungsgemeinschaft under the Contract No. Bel634/1-1, Po256/2-l, and Schw510/2- 
1, respectively. This work was supported by the European Community under contract 
ERBFMGECT950083. 



5 



REFERENCES 



c 



Present address: Fachbereich Physik, Bergische Universitat, D-42119 Wuppertal, Ger- 
many 

^ Present address: Dept. of Physics, Yale University, New Haven CT 06512, USA 

Present address: Max-Planck-Institut fiir Kernphysik, D-69117 Heidelberg, Germany 
Present address: Max-Planck-Institut fiir Physik, D-80805 Miinchen, Germany 
D.H.E. Gross, Rep. Prog. Phys. 53, 605 (1990) 
J. P. Bondorf et al, Phys. Rep. 257, 133 (1995) 

For a recent review of the field see, e.g.. Proceedings of the International Workshop 
XXVII on Gross Properties of Nuclei and Nuclear Excitations, Hirschegg, Austria, 
edited by H. Feldmeier et al. (GSI, Darmstadt, 1999) 
A.S. Botvina and I.N. Mishustin, Phys. Lett. B 294, 23 (1992) 
H.W. Barz et al, Nucl. Phys. A 561, 466 (1993) 
Bao-An Li et al, Phys. Lett. B 303, 225 (1993) 
J. Konopka et al, Prog. Part. Nucl. Phys. 30, 301 (1993) 
Zheng Yu-Ming et al, Phys. Rev. C 53, 1868 (1996) 
A. Schiittauf et al, Nucl. Phys. A 607, 457 (1996) 

Hongfei Xi et a/., Z. Phys. A 359, 397 (1997); Eur. Phys. J. A 1, 235 (1998) 
J. Lauret et al, Phys. Rev. C 57, R1051 (1998) 
S.P. Avdeyev et al, Eur. Phys. J. A 3, 75 (1998) 
Al.H. Raduta and Ad.R. Raduta, Phys. Rev. C 61, 034611 (2000) 
A.S. Botvina et al, Nucl. Phys. A 584, 737 (1995) 
A.S. Hirsch et al, Phys. Rev. C 29, 508 (1984) 
J. Pochodzalla et a/., Phys. Rev. C 35, 1695 (1987) 
D.J. Morrissey et al, Annu. Rev. Nucl. Part. Science 44, 65 (1994) 
for a review of the extensive work done in the intermediate-energy domain see, e.g., 
H. Fuchs and K. Mohring, Rep. Prog. Phys. 57, 231 (1994) 
W. Bauer, Phys. Rev. C 51, 803 (1995) 
K.B. Morley et al, Phys. Rev. C 54, 737 (1996) 
G. Wang et al, Phys. Rev. C 60, 014603 (1999) 
J.A. Hanger et al, Phys. Rev. C 57, 764 (1998) 
K.G. Gulamov et al, Phys. of Atomic Nuclei 62, 533 (1999) 
P.B. Gossiaux et al, Nucl. Phys. A 619, 379 (1997) 
A.S. Goldhaber, Phys. Lett. 53B, 306 (1974) 

for a review of the early pioneering work see J. Hiifner, Phys. Rep. 125, 129 (1985) 
V. Serfling et al, Phys. Rev. Lett. 80, 3928 (1998) 
S. Fritz et al, Phys. Lett. B 461, 315 (1999) 
T. Odeh, PhD thesis, Universitat Frankfurt, 1999 
C. GroB, PhD thesis, Universitat Frankfurt, 1998 
S. Albergo et al, II Nuovo Cimento 89A, 1 (1985) 
J. Pochodzalla et al, Phys. Rev. Lett. 75, 1040 (1995) 
T. Gaitanos et al, preprint |nucl-th/9905020i 1999 



T.C. Sangster et al, Phys. Lett. B 188, 29 (1987) 

B. Berthier et al, Phys. Lett. B 193, 417 (1987) 

E. Renshaw Foxford et al, Phys. Rev. C 54, 749 (1996) 



6 



[37] M. Begemann-Blaich et al, Phys. Rev. C 58, 1639 (1998) 
[38] M. Colonna et al, Phys. Rev. C 55, 1404 (1997) 
T. Lefort et al, preprint |nucl-ex/99100T7| , 1999 



7 



FIGURES 




1 — — r 



1 = 



= * 



o 







u 







o 



I 



i 1 + + ^ • + 



M * ^ * ^ 



• ^\ 

• oj 



30 



o 



20 E 
< 





30 
20 



V 



10 5: 





80 



60 f 
40 ^ 
20 



50 100 

Energy (MeV) 



20 40 

Z 



60 



80 







bound 



FIG. 1. Proton energy spectra measured at Oiab = 150° (left panel) and slope temperatures, 
mean kinetic energies, and relative yields of the two components with high (open circles) and low 
(closed circles) temperatures (right panels). The full lines represent the fit results, the dashed and 
dotted lines give the high and low temperature components for 70 < Z^^^d < 80 individually. 
Note that the fitting was performed with the assumption that two components exist which, for 
Zbound < 30, is not unambiguously supported by the data. 
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FIG. 2. Energy spectra of light charged particles and fragments with Z < A integrated over 
20 < Z^ound ^ 60 {6 lab = 150°). The full lines represent the fit results, the dashed lines indicate 
the trigger threshold of the 300-/Ltm detector. The same normalization is used for all spectra. 
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FIG. 3. Slope temperatures for isotopically resolved charged particles and fragments, measured 
at Qiab = 150° and integrated over 20 < Zhound < 60, as a function of the mass number A. For 
protons and a particles the value of the high-temperature component is shown. 
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FIG. 4. Slope temperatures for light charged particles of mass 2 < A < 4 (squares, triangles 
and diamonds) and the isotope temperature TneLi (dots) as a function of Zhound- The lines give 
the predictions for the fast breakup of a Fermi gas with finite temperature TjjeLi and with densities 
p/po = 1-0 (full line) and 0.3 (dashed). 
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